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Abstract 

Ambartzumian et.al. suggested that the modified Steiner action functional had desirable proper- 
ties for a random surface action. However, Durhuus and Jonsson pointed out that such an action led 
to an ill-defined grand-canonical partition function and suggested that the addition of an area term 
might improve matters. In this paper we investigate this and other related actions numerically for 
dynamically triangulated random surfaces and compare the results with the gaussian plus extrinsic 
curvature actions that have been used previously. 



There has been considerable recent interest in the theory and simulation of random surfaces, inspired 
both by string theory and the study of membranes in solid state physics. The Polyakov partition function 
JlJ for a string embedded in euclidean space with a fixed intrinsic area worldsheet discretizes to 



, N-l 



Z = E/ U dx ^M-s g ), (i) 

where ^ T is a sum over different triangulations of the worldsheet with the same number of nodes N 
which is the discrete analog of the path integration over the intrinsic metric in the continuum. The 
discretized action, S g , is just a simple gaussian 

s 9 = ~ £W-*jT, (2) 

<ij> 

where the variables X{ live on the nodes of the triangulation and the sum < ij > runs over its edges. 
To a solid state physicist this would be the action for a fluid (because of the sum over triangulations) 
membrane with gaussian interactions and no self-avoidance. The difficulties encountered in analytical 
calculations for string theories in physical dimensions are mirrored in simulations of equ. (§) which 
generate very crumpled surfaces from which it is impossible to take a continuum limit jj). Variations on 
the gaussian action, such as using the sum of the area of the triangles making up the surface or the sum 
of their edge lengths, also fail to generate a smooth continuum limit Q],[pj. 

The addition of an extrinsic curvature term, or "stiffness" , to the gaussian action which originally 
arose in the context of biological membranes and QCD || has been mooted as a possible resolution of 
these problems. One discretization of this has shown particular promise, 

S e = (1-ni-nj), (3) 

Ai,A 3 - 

where the normals n^, rij are on adjacent triangles A^, A.,-. The strategy is to examine the phase structure 
of the gaussian plus extrinsic curvature action S g + XS e (henceforward GPEC) as the coupling A is varied 
to see if any continuous transitions are present at which one might define a non-trivial continuum theory 
The case of fixed, or "tethered" , surfaces where there is no sum over triangulations offers some hope 
because there does appear to be a second order transition between a low A crumpled phase and a large 
A smooth phase [] for these when there is no self-avoidance J ^J. Initial work on dynamical surfaces with 
modestly sized meshes tended to support a similar picture |9|, with a peak in the specific heat that grew 
with mesh size. More recent simulations [jn^jllj with higher statistics and larger meshes have cast doubt 
on this because the increase in the peak height tails off with increasing mesh size. This is not in itself 
worrying since it is possible that the string tension and mass gap may still scale in such a way as to give 
a non trivial continuum limit, even at a higher order transition JTo| | . However, it was pointed out in Jll| 
that the data might also be construed as indicating a rapid crossover rather than a true transition or 
even the presence of a low mass bound state like that in the 2d 0(3) model. 

This rather confused picture for GPEC actions naturally prompts the question of whether discrete 
models with more clearly cut phase transitions exist. It was suggested in jl^] that an action based on 
the modified Steiner functional [^3| might be a suitable alternative candidate as it possessed the nice 
geometrical property of being subdivision invariant as well as retaining a desirable stiffening effect on the 
surfaces. The action was 

w = ji;w-^i9w 1 (4) 

<ij> 

where O(cnj) — |tt — | and onj is the angle between the embedded neighbouring triangles with common 
link < ij >. This is essentially a coarse discretization of the absolute value of the trace of the second 
fundamental form of the surface. However it was observed in E3 that an action of this form ran into 
problems with the entropy of vertices in smooth configurations and failed to give a well-defined grand 
canonical partition function. Our previous microcanonical simulation of such an action p"5[ ] had shown 
that it produced smooth surfaces but in this case of course it is not clear how to obtain the continuum 
limit, where the problems of |l4f| would presumably resurface. 

1 These have been proposed as models of polymerized membranes in solid state physics j§j 
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We had observed in fig] that varying the coupling A in an action of the form 

si^E \ x " - x i\ + A E \ x * ( 5 ) 

<ij> <ij> 

would allow one to employ the approach used in our earlier simulations of GPEC actions, namely hunting 
for continuous transitions at which to define the continuum theory. We chose the first, edge-length, term 
in S± because of its simplicity and similarity to the second Steiner term. The authors of |l4| pointed out 
that an subdivision invariant action of the form 

S 2 = £|A|+A£ W-XfMav), (6) 

A <ij> 

where | A| is the area of a triangle, might cure the entropy problems of S s t e iner by virtue of the first term 
being coercive enough to compete with the entropy of smooth surfaces. Si might also be expected to 
effect a cure for the same reason, but it is no longer subdivision invariant. It is not clear that this is 
a prerequisite for a well-behaved action in any case as the GPEC action does not share this property 
either. If we do not insist on subdivision invariance we could also envisage combining a gaussian term 

with S s teiner 

S3 = \ E x rf + A E \ x i X WM- (7) 

<ij> <ij> 

The arguments presented in |l2] ] suggested that S s t e iner was essentially a discretization of the absolute 
value of the trace of the second fundamental form, at least for surfaces embedded in three dimensions. 
Making use of the original definition of the trace of the second fundamental form 

K = Kl~ d a X» ■ d a n^ (8) 

we could also attempt a direct discretization in the three dimensional case. Various possibilities are 
indicated in Fig.l. One could take 

\K\= J2 |(*i-*i)-(ni-n2,3)|. (9) 

where the subscript 2, 3 denotes a choice of either normal 712 or 113 in the discretization 0of<9n. If we 
pick ni it disappears in the sum because it is perpendicular to Xi — Xj , whereas a coarser discretization 
of the derivative given by choosing retains both of the normals ni and 77,3. We can then pair these 
alternative discretizations with either edge-length, area or gaussian terms to give the further possibilities 



s * = \ E \xt-x^\+xj2 \ R \ 

<ij> <ij> 

5 5 = ]T|A| + A]T \K\ 

A <ij> 

5 6 = ^(Xi-Xjf + X^IK]. (10) 

<ij> <ij> 

In this paper we carry out a qualitative exploration of the phase structure of Si , £2 , S3 , S4, S5 and Sq 
to see if there is any justification for expecting a well-behaved theory. In order to do this we measured 
the standard repertoire of intrinsic and extrinsic variables. We included a local factor in the measure for 
compatibility with our earlier simulations which can be exponentiated to give 

S m = ^logfe), (11) 

i 

where qi is the number of neighbours of point i, and d = 3 dimensions. If we denote the general form 
of actions Si, ...,S , 6 by S a + \Sp we thus simulated S a + XSp + S m . We measured < S m > and the 

2 Choosing Xj — instead of Xj — Xj gives identical results 
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mean maximum number of neighbours < max(qi) > to get some idea of the behaviour of the intrinsic 
geometry. The extrinsic geometry was observed by measuring < Sp > and its associated specific heat 



a 



= ^(<sl>-<s >' 



(12) 



as well as the gyration radius X2, a measure of the mean size of the surface as seen in the embedding 
space, 

X2 = V (X? - X?) 2 q iqi . (13) 

9N(N - 1) jr 1 V 1 0> 3 y ' 

Note that the sum for X2 is over all pairs of X's on the mesh. The expectation value of the first term 
in the action < S a > was also measured for completeness, although unlike the GPEC actions we can no 
longer use its value to confirm that equilibration has occurred. 

The simulation used a Monte Carlo procedure which we have described in some detail elsewhere Jl6| . 
It first goes through the mesh moving the X's, carrying out a Metropolis accept/reject at each step, and 
then goes through the mesh again carrying out the "flip" moves on the links, again applying a Metropolis 
accept/reject at each stage. The entire procedure constitutes a sweep. Due to the correlated nature of 
the data, a measurement was taken every tenth sweep and binning techniques were used to analyse the 
errors. We carried out 10K thermalization sweeps followed by 30K or 50K measurement sweeps for each 
data point. The acceptance for the X move was monitored and the size of the shift was adjusted to 
maintain an acceptance of around 50 percent. The acceptance for the flip move was also measured, but 
in this case there is nothing to adjust, so as for GPEC actions this dropped with increasing A (but was 
still appreciable even for quite large A). 

We now move on to consider the numerical results in more detail for action Si (edge-length plus 
Steiner term) on a 72-node surface. We can see from Table. 1 
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Table 1 
Results for S t , N=72 

that there does, indeed, appear to be a transition to a smooth phase at large A, as inspection of Figs. 2, 3 
confirms. We can see in Fig. 2 that there is a peak in the associated specific heat Cp at around A = 3.25 
where we also observe a sharp jump in the measured size of the surface as given by X2 which is plotted 
in Fig. 3. The smooth nature of the large A phase is corroborated by the snapshot of a surface at A = 5 in 
Fig. 4. For small A the surfaces are crumpled as can be seen in Fig. 5. The behaviour in the crumpled phase 
is, however, different from that seen with GPEC actions. Firstly, at A = the surfaces are larger than 
those seen with a purely gaussian action, presumably because the edge length action is less confining than 
the gaussian action. Secondly, and rather strikingly, X2 decreases with increasing A up to the putative 
phase transition point, where it increases again. This can be understood by noting that the Steiner term, 
Sp, is not scale invariant like the extrinsic curvature. There are thus two ways of decreasing the value 
of Sp, one can either make the surface smaller to decrease the |JQ — Xj\ terms, or make the surface 
smoother to decrease the 9 terms. It is evident that the first occurs for A up to the phase transition, as 
can be confirmed by looking at S a , which is a measure of the sum of \Xi — Xj \ terms - this decreases. The 
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behaviour of X2 in the crumpled phase appears to be a generic feature of all of the actions containing 
Steiner terms as it is also observed for S2 (area + Steiner) and S3 (gaussian + Steiner) actions. The 
Steiner term Sp appears to correlate well with the extrinsic curvature, as can be seen in Table. 1, where its 
value is high in the crumpled low A phase (where the extrinsic curvature, which we have not tabulated, is 
also high) and drops off as A increases and the surfaces become smoother. The behaviour of the intrinsic 
observables S m and max(qi) is very similar to that of the GPEC actions: S m increases smoothly through 
the phase transition and max(g^) decreases, indicating that the intrinsic geometry becomes slightly more 
regular at larger A. 

52 (area plus Steiner term) gives the results in Table. 2. Again we have a smooth phase at large A 
with a transition to a crumpled phase at lower A. Snapshots of surfaces in the smooth phase look similar 
to that produced by S\ in Fig. 4 and the low A phase is, if anything, even spikier than that of Si. In 
fact, at A = the surfaces appear to be unstable to the formation of very large spikes, appearing as a 
collection of thin whiskers emanating from a central point. We have not included the values for A = in 
Table. 2 and Figs. 2, 3 as X2 and Sp are very large. As we can see from Table. 2 
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Table 2 
Results for S 2 , N = 72 

the X2 values for S2 even with non-zero A are larger for a given coupling than for 5*1, suggesting the 
spikiness persists away from A = 0. We can also see from the table that X2 decreases at small A as does 
S a for precisely the reasons indicated above for £1 . It is clear that the transition is taking place at a lower 
A value than for Si, as both the peak in the specific heat and the jump in X2 are at A ~ 2.25, rather 
than A ~ 3.25. At first sight this might appear rather surprising as we have just seen that the crumpled 
phase is spikier for S2 rather than Si, so we might expect to have to work harder to escape from it and 
thus have a phase transition at larger A. However, the dimensions of the area and edge length terms are 
different so we cannot directly attach significance to the numerical value of the Steiner coupling. We can 
see in Fig.l that the peak in the specific heat for S2 appears to be considerably stronger than for Si []. 

S3 (gaussian plus Steiner term) gives the results shown in Table. 3. where we can see that the behaviour 
is similar to S2. Although S3 is not subdivision invariant there is still a larger peak in the specific heat 
than for Si as can be seen clearly in Fig. 2. The peak is at a larger A (~ 3) value than for S2, which 
is perhaps surprising as one might have expected that it would be easier to escape from the crumpled 
phase of S3 because it is less "spiky" than that of S2. On the other hand, as can be seen in Fig. 3, the X2 
values of S3 are more similar to Si than S2 . The characteristic decrease in X2 for small A is also present 
for S3. Again the values for A = are omitted, this time because they are identical to the GPEC action 
at A = 0. 

3 It is worth remarking that the sharpness of the peak in the specific heat should not be regarded as the criterion for a 
good candidate continuum theory as models which do have second order transitions on fixed lattices, such as Ising or Potts 
models, have higher order transitions on dynamical lattices with well behaved continuum limits. The GPEC action j^and 
Si may be similar. 
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Table 3 
Results S 3 ,N = 72 



We might naively expect S4, S5, Sq to behave in a similar fashion to S\, S2, S3 respectively, but a glance 
at Fig. 6 and Fig. 7 where we show snapshots of 72-node surfaces generated by S4 and S5 respectively at 
A = 5 reveal that this is not the case. It is clear for both of these that Sp is failing to smooth out the 
surfaces at large A, for S4 (and Sq too, which we have not shown) because of large numbers of "back to 
back" triangles and for S5 because it fails to prevent the surface degenerating into thin whiskers, similar 
to those observed for S2 at A = 0. In Table. 4 
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Table 4 
Results for £4, N=36 



the numerical results for S4 show that despite the obvious roughness of the surfaces at large A, Sp which is 
supposed to discretize the absolute value of the trace of the second fundamental form, is decreasing with 
increasing A. The numerical results for S5, Sq are similar to this - although Sp decreases with increasing 
A the surfaces obviously remain rough. We have not pursued simulations of S4, S5, Sg on larger surfaces 
because of the obvious pathologies seen in these results. 

For the choice of in the discretization of \K\ in equ.@ the nature of the pathology can be pinned 
down by looking at Sp without the modulus sign for the various surfaces. Summing up the various 
contributions on the smooth surfaces generated by Si, S%, S3 gives a much larger total than summing up 
the contributions on surfaces generated by S4 , S5 , Sq where we see a large degree of cancellation between 
terms with alternating signs. The reason for this is clear if we examine equ.(|9|): there is an ambiguity in 
the direction of rii if it is forced to be perpendicular to X,- t — Xj. It can be at either +90° or —90° and the 
Monte Carlo will pick both with equal probability, leading to normals reversing direction from triangle 
to triangle and the non-smooth surfaces that are observed with S4, S5, Sq. This problem does not arise 
for Si, S2, S3 where the Steiner term is minimized when the angle between adjacent triangles is 180° and 
there is no possibility of confusion in the orientation of the normals. Choosing 713 in the discretization of 
I if I in equ.(||), which retains both normals, might be expected to give better results as it is more difficult 
to generate crumpled configurations with both ni and 713 perpendicular to Xi — Xj. Nonetheless, this 
alternate discretization also fails to produce a smooth phase. Examining Sp without the modulus sign 
again reveals a large degree of cancellation due to alternating signs in the surfaces generated by Si, S5, Sq 
(which is not seen for the smooth surfaces produced by Si, S2, S3). 

The simulations described in this paper can be regarded as a qualitative study of the feasibility of 
using alternatives to the GPEC action for random surfaces. We have found that Si and S3 which are 
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not subdivision invariant and S2, which is, all seem to be good candidates for further exploration. S4, 
S5 and 56 on the other hand, fail to discretize the trace of the second fundamental form properly. It 
thus appears that the choice of "stiffness" term essentially determines at least the qualitative behaviour 
of the surface at the crumpling transition. Similar conclusions have been drawn by one of the authors 
for various versions of the extrinsic curvature combined with area, gaussian or tethering terms on fixed 
random surfaces §L7\ . The crumpled phase for the various actions, particularly for 5*2 at very small A, 
does display variations. The interesting question of whether the subdivision invariance of S2 has an effect 
on the random surface model as it approaches the continuum limit has not been investigated in this 
paper because we have made no attempt to carry out finite size scaling for Si, 52, 53 at the transition. 
This would also allow a more quantitative comparison with the GPEC actions and establish whether the 
properties of the transition for Si,S2, 53 are more clear cut than those seen in [^0), p| . These issues, 
and the equally important question of whether the string tension and mass gap scale, are currently being 
addressed in a much larger simulation. 

The comparison with the results of the GPEC action will also shed some light on the still murky 
problem of universality for random surfaces |Q. As we have seen, there is a wide choice of possible 
terms and discretizations for a random surface action and the question of whether they lead to the same 
continuum theory has not been clearly answered. In the longer run the best approach may be to use 
the Monte Carlo renormalization group, perhaps starting on rigid meshes because of the computational 
complexity of dynamical meshes. This would then allow consideration of all possible terms up to a given 
weight. We are currently exploring the feasibility of this. 
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and an Spanish/British Acciones Integradas grant (DE and DAJ). CFB is supported by DOE under con- 
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Figure Captions 

Fig. 1. The edges and normals that are involved in the discretization of the the Sp term in S4, £5 and 
S e . 

Fig. 2. The specific heat Cp for actions Si, 52, S3. 
Fig. 3. The gyration radii X2 for actions Si, S2, S3. 
Fig. 4. A snapshot of a mesh generated by Si with A = 5. 
Fig. 5. A snapshot of a mesh generated by Si with A = 0. 
Fig. 6. A snapshot of a mesh generated by S4 with A = 5. 
Fig. 7. A snapshot of a mesh generated by S5 with A = 5. 
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